Toward a Midisuperspace Quantization of LeMaitre-Tolman-Bondi Collapse Models 
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LeMaitre-Tolman-Bondi models of spherical dust collapse have been used and continue to be used 
extensively to study various stellar collapse scenarios. It is by now well-known that these models 
lead to the formation of black holes and naked singularities from regular initial data. The final out- 
come of the collapse, particularly in the event of naked singularity formation, depends very heavily 
on quantum effects during the final stages. These quantum effects cannot generally be treated 
semi-classically as quantum fluctuations of the gravitational field are expected to dominate before 
the final state is reached. We present a canonical reduction of LeMaitre-Tolman-Bondi space-times 
describing the marginally bound collapse of inhomogeneous dust, in which the physical radius, R, 
the proper time of the collapsing dust, r, and the mass function, F, are the canonical coordinates, 
R(r), r(r) and F(r) on the phase space. Dirac's constraint quantization leads to a simple functional 
(Wheeler-DeWitt) equation. The equation is solved and the solution can be employed to study some 
of the effects of quantum gravity during gravitational collapse with different initial conditions. 
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I. INTRODUCTION 

In the standard treatment of Hawking radiation 
from black holes, one begins with a black hole that is 
formed in some classical model of gravitational collapse, 
surrounds the hole by a quantum field and examines the 
propagation of this field on the classical background pro- 
vided by the hole. The quantum field behaves as a ther- 
mometer. The quantum modes within the horizon are 
averaged over and one finds that the black hole radiates 
thermally leading to its "evaporation". The tempera- 
ture that characterizes the evaporation of a black hole 
is inversely proportional to its mass and the radiation 
flux is inversely proportional to its mass square. For an 
astrophysical black hole, the radiation flux and temper- 
ature are therefore so small that the semi-classical ap- 
proximation is expected to be an adequate description of 
the evaporation until the hole is roughly of Planck di- 
mensions at which point higher order quantum gravity 
effects will undoubtedly become important. This means 
that the final state of the black hole will depend on quan- 
tum gravity. The black hole may evaporate completely 
or it may leave a remnant. If it evaporates completely, 
it is important to understand what happens to the infor- 
mation that was initially trapped within its horizon. 

Classical models of collapse also lead to the formation 
of naked singularities for regular initial data 0. One 
may then ask if the quantum radiation from naked sin- 
gularities is similar to the radiation from black holes. By 
surrounding a classical naked singularity with a quan- 
tum field and examining its quantum modes at null in- 
finity one finds that the evaporation of a naked singu- 
larity is qualitatively distinct from that of a black hole 



The radiation flux diverges as the Cauchy hori- 
zon j7|-pT| is approached and the spectrum of the radi- 
ation is non-thermal [fslJl^l, falling off as the inverse fre- 
quency. Contrary to the case of an astrophysical black 
hole, the (divergent) flux should be observable and the 
unique spectrum should serve to distinguish objects un- 
dergoing this type of collapse from other celestial emit- 
ters. However, a closer look at the semi-classical approx- 
imation just described reveals that the flux of radiation 
is essentially negligible (on the order of one Planck mass) 
until about one Planck time before the putative Cauchy 
horizon is reached |l| . The semi-classical approximation 
therefore signals a quantum instability of naked singular- 
ities (and therefore a mechanism for the Cosmic Censor) 
but its quantitative predictions must be tested in a full 
quantization of all the degrees of freedom, including the 
gravitational field. When quantum gravitational effects 
are accounted for, does the flux continue to diverge and 
the radiation spectrum continue unique? Indeed, does 
the quantum theory serve as a Cosmic Censor and are 
there significant observational consequences of collapse 
into naked singularities as the semi-classical approxima- 
tion suggests? 

In a step toward answering these and other ques- 
tions regarding the final stages of gravitational collapse 
Kuchaf |L4j examined a midi-superspace quantization 
of the Schwarzschild black hole, presenting in the pro- 
cess a remarkable series of canonical transformations 
that greatly simplified the dynamical equations. In the 
present paper, we describe a generalization of this work 
to spherically symmetric, marginally bound LeMaitre- 
Tolman-Bondi space-times. We show below that there 
is an analogous description of the gravitational part of 



the action in terms of the "mass function" , the physi- 
cal radius and their conjugate momenta. Furthermore, 
the coupling to dust introduces the proper time of the 
collapsing matter and its conjugate momentum. Thus 
time evolution appears naturally into the dynamical con- 
straints. 

The hypersurface action yields two constraints, viz., 
the Hamiltonian constraint and the momentum con- 
straint, which are given in terms of this canonical chart 
consisting of the mass, F[p(r)], contained within spheri- 
cal shells of fixed shell-label coordinate, p{r), the physical 
radius, R(r), the dust proper time, r(r), and their con- 
jugate momenta, Pp(r), Pr{t) and P T (r) respectively. 
Here r is the radial label coordinate of a foliation of the 
space-time by spacelike hypersurfaces. The momentum 
conjugate to the mass function, Pp(r), may be elimi- 
nated in the Hamiltonian constraint using the momentum 
constraint. This leads to a new and simpler constraint 
that is able to take the place of the original Hamiltonian 
constraint. Dirac's constraint quantization then yields 
a two dimensional Klcin-Gordon-likc functional equation 
with a potential term that depends on the mass function 
and its derivative with respect to the label coordinate r. 
The simplest possible scenario, in which the mass func- 
tion is constant throughout the space-time, the same for 
all shells, describes the Schwarzschild black hole. The 
Schwarzschild black hole will thus emerge as a special 
case of the general class of models we quantize below. 

In section II we review the classical collapse models 
being considered in this article and present the canoni- 
cal formulation for spherically symmetric space-times in 
section III, where we also discuss the fall-off conditions 
appropriate to the models and the resulting boundary 
terms. We reconstruct the mass and time from the canon- 
ical data in section IV. This leads naturally to new vari- 
ables viz., the mass, the dust proper time, the physical ra- 
dius and their conjugate momenta which are introduced 
along with the generator of the canonical transformation 
from the old to the new variables. We apply Dirac's 
quantization program to the new constraints in section 
V, showing how the new variables lead to a simplified 
Wheeler-DeWitt equation. We then present a solution 
for an arbitrary, but non-constant mass function. 



II. CLASSICAL MODELS 

The LeMaitre-Tolman-Bondi models ]L5| constitute a 
complete solution of the Einstein equations for a matter 
continuum of inhomogeneous dust, i.e., they are solutions 
of the spherically symmetric Einstein's field equations, 
G M „ = with vanishing cosmological constant 

and with stress-energy describing inhomogeneous, pres- 
sureless dust given by — eu^u v . The solution is the 
LTB metric, given in co-moving coordinates as ]16| 



ds 2 = dr 2 - 



R 2 



1 + / 



dp 2 - R 2 dtt 2 , 



F / F 



(2.1) 



(we have set 8ttG = 1 = 6) where R is the physical ra- 
dius. A tilde (~) represents a derivative with respect to 



p and a star (*) represents a derivative with respect to 
the dust proper time r. The functions f(p) and F(p) are 
arbitrary functions only of p, interpreted respectively as 
the energy and mass functions. The energy density of the 
collapsing matter is e(r, p), and the negative sign in the 
third equation above is required to describe a collapsing 
cloud. Its general solution is given up to an arbitrary 
function ip(p) of the shell label coordinate. This arbi- 
trariness reflects only a freedom in our choice of units 
i.e., at any given time, say r Q , the function R(r ,p) can 
be chosen to be an arbitrary function of p. 

The mass function, F(p), represents the weighted mass 
(weighted by the factor + /) contained within the 
matter shell labeled by p. If a scaling is chosen so that the 
physical radius coincides with the shell label coordinate, 
p, at t = 0, then it can be expressed in terms of the 
energy density at r = according to 



Hp) 



e{0,p)p 2 dp, 



(2.2) 



while the energy function, f(p), can be expressed in terms 
of the initial velocity profile, v(p) = R*(0,p), according 
to 

f{p) = v 2 (p) --J e(0,p)p 2 dp. (2.3) 

The marginally bound models, which we will consider 
in this paper, are defined by f(p) — 0. For the scaling 

referred to a bov e, we must choose ip(p) = p^ , whence the 
solution of (2.1) can be written as 



RHt,p) 



(2.4) 



The epoch R = describes a physical singularity, whose 
singularity curve 



r(p) 



2pi 



(2.5) 



gives the proper time when successive shells meet the 
central physical singularity. Various models are obtained 
from choices of the mass function, F(p). For example, the 
Schwarzschild black hole is the marginally bound solution 
with F(p) = 2M, a constant. 

A collapsing star does not have sharp boundaries. In 
the simplest possible approximation to reality, it will con- 
sist of a dense core surrounded by a crust of lower den- 
sity which, itself, is encased in a cloud whose density will 
smoothly go to zero with distance from the star's cen- 
ter |17]]. Nevertheless, as an approximation, one could 
consicTer a sharp boundary at some constant shell label, 
Pb, with an exterior Schwarzschild metric. One could 
also consider several regions described by different mass 
functions describing successively lower matter densities 
as one moves out from the center. The mass function as 
a whole would not be differentiable at the boundary or 
boundaries, but we require it to be continuous. 

One can then examine outgoing families of non- 
spacelike geodesies and check if there exist congruences 
that terminate in the past at the central singularity [fl8| . 
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If such congruences exist then the collapse leads to a 
naked singularity and if they do not exist then the col- 
lapse leads to a black hole. A naked singularity may fur- 
ther be characterized as globally naked or locally naked 
depending on whether the outgoing geodesies succeed in 
reaching null infinity or not. 

In general one finds that both black holes and naked 
singularities may develop as the end states of collapse, 
depending on initial data, i.e., on the initial density and 
velocity profiles of the collapsing dust. For example, in 
the marginally bound self-similar collapse model, F(p) — 
A/9, where A = const., both outcomes described by the 
Penrose diagrams in figure 1 are possible, depending on 
whether A > 0.1809 (black hole) or A < 0.1809 (naked 
singularity) . 



" Event Horizon 
Cauchy Horizon 




Fig.l: Black hole (left) or naked singularity (right) 

It is in fact believed by many that sets of initial data may 
evolve in general relativity toward either naked singular- 
ities or black holes independently of the equation of state 
or the type of matter used. 

Geometrodynamics views space-time as the dynami- 
cal evolution of spatial hypersurfaces. When the collapse 
evolves toward a black hole such spatial hypersurfaces ex- 
ist, starting at infinity, crossing the horizon and contin- 
uing to p — without encountering the central singular- 
ity. On the contrary, when the collapse evolves toward a 
naked singularity, the spatial hypersurfaces in the future 
of the initial singularity cross the Cauchy horizon and col- 
lide with the central singularity. No sensible boundary 
conditions can be specified on a singularity and evolution 
in the future of the initial singularity is arbitrary. To 
avoid the consequent breakdown of predictability, Pen- 
rose proposed the Cosmic Censor |L9| which, as men- 
tioned in the introduction, is likely to be the quantum 
theory itself and will come into play before the Cauchy 
horizon has a chance to form. It is of particular inter- 
est, therefore, to understand precisely how the system 
behaves close to, but in the past of, the putative Cauchy 
horizon. Spatial hypersurfaces in the past of the Cauchy 
horizon are well defined and the quantum evolution of the 
system may be studied until the time of its formation. 



ized by two functions, L(r) and R(r) of the radial label 
coordinate, r, according to 



da 2 = L 2 (r)dr 



R 2 (r)dQ 2 



(3.1) 



where f2 is the solid angle. Neither L(r) nor i?(r) can be 
negative and we take them to be positive definite except 
possibly at the center. i?(r) represents the physical ra- 
dius of the point labeled by r on the surface. It behaves as 
a scalar under transformations of r, whereas L(r) behaves 
as a scalar density. The corresponding four dimensional 
line element may be written in terms of two additional 
functions, the lapse, N(t,r), and the shift, N r (t,r), as 



ds 2 



N 2 dt 2 - L 2 (dr- N r dt) 2 - R 2 dn 2 . (3.2) 



In this spherically symmetric space-time, we will consider 
the Einstein-Dust system described by the action 



S = 



-L / d 4 xV=g n 

l07T 



i- J d 4 x^ e(x) [g a pU a U? + l] 



(3.3) 



where 1Z is the scalar curvature. As is well known, the 
gravitational part of this action can be cast into the form 



S 9 = dt 



dr 



P L L + P R R - NH 9 - N r H? 



C9 



(3.4) 



with the momenta conjugate to L and R respectively 
given by 



Pl 
Pr 



R 

N 
1 

N 



-R + N r R' 
-LR- LR+ (N r LR)' 



(3.5) 



and where the overdot and the prime refer respectively to 
partial derivatives with respect to the label time, t, and 
coordinate, r. The lapse, shift and phase-space variables 
are required to be continuous functions of the label co- 
ordinates. The boundary action, Sg S , is required to can- 
cel unwanted boundary terms in the hypersurface action, 
thereby ensuring that the hypersurface evolution is not 
frozen. It is determined after fall-off conditions appropri- 
ate to the models under consideration are specified. The 
super-Hamiltonian and super-momentum constraints are 
given by 



H 9 



Pt.P. 



LP 2 



m = r'p r 



r 

L R' 
2 ~ ~2L 

LP'l 



2R 2 

i 



( RR' 



(3.6) 



III. CANONICAL DYNAMICS 

The line element da on a spherically symmetric three- 
dimensional Riemann surface £ is completely character- 



We will assume that the chosen mass function, F(o), 
is such that at infinity Kuchaf's fall-off conditions Ej] 
are suitable and we will adopt them here. These con- 
ditions would be applicable, for example, in models in 
which the collapsing metric asymptotically approaches or 
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is smoothly matched to an exterior Schwarzschild back- 
ground at some boundary, p^. They read 



L{t,r 
R(t,r 

Pn(t,r 
N(t,r 
N r (t,r 



1 - 
r - 
0°°(r 
0°°(r 
N+(t) 
0°°(r 



M + {t)r~ l 

o°°( r - e ) 



(3.7) 



Again, because the label radial coordinate r € [0, oo), we 
must also consider the boundary conditions at r = 0. Let 
the mass function near the center (p = 0) have a series 
expansion of the form 



F(p) = F + F lP + F 2 p 2 + ... = Y,FnP r< 



(3.8) 



If F n = V n > but F a > 0, the solution describes the 
Schwarzschild black hole of mass F a /2. The marginally 
bound, self-similar model mentioned in the previous sec- 
tion corresponds to F\ = A > and F n = V n ^ 1, 
or a density profile that behaves as e(0, p) ~ p -2 . Thus 
t = is the singular epoch for the self-similar model and 
this singular density profile arises from a regular initial 
profile at some r < 0. 

Consider models with F(p) = F n p n . As far as the fall- 
off conditions at the center are concerned, two classes of 
models arise, viz., n < 3 and n > 3. The considerations 
below are applicable to models with n < 3, although 
conditions appropriate to models wi th n > 3 may likewise 
be given. Referring to equation (2.4) we find that as 
p ; R{t, P ) r> i? (r)pf + G(pt+ fe ) and L(t, p) « 
L (r)pt- 1 + 0(pi ~ 1+fc ), where k > 0. We will exclude 
the black hole (n = 0) in the following because in that 
case the space-time may be analytically continu ed t o r — 
— oo where the boundary conditions given in (3.7) may 
be applied and no conditions at r = need be given. For 
a genuine collapse, the cases of interest are those with 
ri > around the central region. 

Let us assume that, at the center, (p, r) approach (r, t) 
as p = r + 0(r a ) {a > 1) and r = f(t)+0(r a ). Now, in 
the following section we will show that the mass function 
is recovered locally from the canonical data according to 



F = R 



Pl_R^ 

1 + R 2 L 2 



(3.9) 



Again, with p = r + 0(r a ) and r = /(*) + 0(r a ), it 
follows that 



R' 



l + O^"- 1 ). 



(3.10) 



and therefore 



F = F n p n = F n r n + 0(r n+a - L ) 

p2 

= -k + 0( r «/3+«-l) ) 

R 



(3.11) 



F„ = -^r' 4 ™/ 3 + C>( r -2n/3+a-l^ (3 12 ) 

R 



The last equation can be satisfied if Pl falls off faster 
than r 2 ™/ 3 and a = 1 + 2n/3. Furthermore, requiring 
both terms in the Liouville form to have the same be- 
havior at the origin, we choose the following conditions 
near r = when n = {1, 2, 3}: 



Ro(t) 



»«/3 



R(r,t) 

L(r,t) = L {t)r n '*- 1 
P L (r,t) = o(r 2n/3+1+e ) 
P R {r,t) =0(r 2n / 3+e ) 
N(r,t) =N (t)+0(r e ) 
N r (r,t) = 0(r e ) 



,n/3+e 



+ 0{r n/3 - 1+e ) 



(3.13) 



The conditions (3.7) and (3.13) ensure that the Liouville 
form is well behaved both at the origin (0(r'™ +e ))) as 
well as at infinity (O 00 (r~ e )) . The Hamiltonian and mo- 
mentum densities behave as 

H 9 = 0{r n+1+e ), Hf, = 0(r"~ 1+e ) (3.14) 

as the origin is approached, and 

H 9 =0°°(r- 2(1+e) ), H 9 = C°°(r- (1+e) ) (3.15) 

asymptotically. Thus the total Hamiltonian and momen- 
tum are well defined and the surface action is meaningful. 
The potential contributions to the s urfa ce action can be 
read off the constraint equations, ( |3.6| ), by considering 
variations of the phase space variables. Applying the 
fall-off conditions at infinity, one finds that only one of 
these is non- vanishing and behaves as 



dtN + (t)6M+(t). 



(3.16) 



On the other hand, with the fall-off conditions in (3.13) as 
r — ► 0, all the variations vanish at the origin. Therefore 
the only contribution to the boundary action comes from 
the surface term at infinity and we find 



dtN+(t)M+(t). 



(3.17) 



We will return to this surface action shortly. 

Let us now turn to the dust portion of the action 
inQ, 



S d = I d^xj~g e(x) [g a pU a U fi + l] . (3.18) 

It has been exhaustively analyzed by Kuchaf and Torre 
p(| and by Brown and Kuchaf pl| . It may be under- 
stood in two ways: either as a consequence of imposing 
coordinate conditions (the Gaussian conditions) or as a 
realistic material medium. For the collapse problem it 
is a realistic material medium and for the LTB models 
being considered, it is non-rotating. Dust is described by 
eight space-time scalars, e, t, Z k and Wk (k € {1, 2, 3}). 
The physical interpretation of these variables which fol- 
lows from an analysis of the equations of motion was 
given in ref. pl[ and will be summarized here for com- 
pleteness, t is the proper time measured along particle 
flow lines, Z k are the comoving coordinates of the dust, 
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W k are the spatial components of the four velocity in 
the dust frame, and e is the dust proper energy den- 
sity. All these scalars are assumed to be functions of the 
space-time coordinates. In particular, the four variables, 
Z K = (r, Z k ), are independent functions, det|Z K M | =/= 0, 
and the four- velocity of the dust particles may be defined 
by its decomposition in the co-basis Z K by 



w k z k 



(3.19) 



In the spherically symmetric geometry described by (1.5), 
the dust action may be cast into the form 



5'" 



where 



dt / dr 



P t t + P k Z k - NH - N r H c r 



Pk 



LW 

-w k p T 



e(r,t) [-U t + N r U r 



(3.20) 



(3.21) 



are the momenta conjugate to the dust proper time and 
the frame variables respectively, and 



H d = P T 



U 2 
1 + — 

L 2 



m = -u r p T . 



(3.22) 



This expression for H d is obtained by exploiting the 
fact that e(t, r) is a Lagrange multiplier, a nd th erefore 
S£/5e = 0. When the dust is non-rotating (3.22) is fur- 
ther simplified by requiring that the dust motion be de- 
scribed with respect to the frame orthogonal foliation. 
Then we may impose the additional constraints P k = 
(when imposed on the state functional these constraints 
simply mean that the state functional does not depend 
on the frame variables Z k ). Using this in ( 3.21 ) we see 
that W k = and so U r — — r'. 

Combining the results thus far obtained, we now give 
t he f ull form of the canonical description of the action in 
(3.3) as: 



S= dt 



dr [P t t + P l L 

+ P R R - NH - N r H r 



•S| E (3.23) 



with 



H= - - 




and 



//,. — 



5 as = / dtN+{t)M+{t). 



(3.24) 



(3.25) 



for the boundary contribution. 



IV. NEW VARIABLES 



The hypersurfaces we consider, from which (3.2) is con- 
structed, must eventually be emb edde d in a space-time 
described by the metric given in (2.1) with / = 0. We 
imagine that they are leaves of the foliation r(t,r) and 
p(t , r). Then the functions L(t, r) and R(t, r) appearing 



in (3.2)are easily determined by substituting the foliation 
in (O). We find 



L 2 



N r 
N 



52 /2 ,2 
Hp — T 

R 2 pp' - tt' 



R 



L 2 

(rp' - pt 1 ) 



(4.1) 



The last of the equations above involves taking a square 
root. We must check that the positive square root taken 
leads to a positive lapse function in all the regions of 
the space-time. Call T = 1 — F/R where F is the mass 



function and note that R* = — yl — T according to (2.1) 
with / = 0. 

Substit utin g the expression for N(t, r) and N r (t, r) 6b- 
) into the expressions for the canonical mo- 
), one obtains the relation 



tained in ( 14.1 
menta in 

LPl 



R 



R(t P ' - T'p) 



-R(R 2 p' 2 -t' 2 ) 

+ R'(R 2 pp' -rr% (4.2) 



which, after some algebra, can be used to obtain r' in 
terms of the canonical variables, 



LP L R'Vl^T 



RT 



(4.3) 



Inserting this into the expression for L 2 in (4.1), we find 
that 



L = 



R' 2 



L 2 Pj 
R 2 T ' 



which determines T , 
T = 



IT 

T 2 R 2 ' 
and, through J 7 , recovers the mass function 

,/2~ 



F = R 



pl ir_ 

+ R 2 L 2 



(4.4) 



(4.5) 



(4.6) 



This is the relation used in the previous section when dis- 
cussing the fall-off conditions at the center. It enables us 
to determine the mass function locally from the canonical 
data. Furthermore, once the dust proper time is fixed at 
some poi nt on the hypersurface, say at spatial infinity, 
equation ( |4.3| ), which determines the difference in dust 
proper times between any two points ri and on a spa- 
tial hypersurface, will determine it at any point on the 
hypersurface. Note that at the horizon, when T = 0, 
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R!_ 
L 



Pl 
R ' 



t' continues to be well behaved, as expected. 
It turns out that the functions Pp, defined byQ 



Pf = 



LPl 
2HF' 



(4.7) 



(4.8) 



and the mass function, P, form a conjugate pair of vari- 
ables. Moreover, because neither F nor Pp depend on 
Pr, they have vanishing Poisson brackets with R. They 
also have vanishing Poisson brackets with r and P T . 
Their Poisson brackets with Pr, however, do not van- 
ish and one cannot directly replace the pair (L, Pi) with 
the more transparent variables (F, Pp) to form a new 
chart. Instead one asks if it is possible to determine 
a new momentum, _Pr, conjugate to R and such that 
the set (r, P, P, P t ,Pr, Pf) forms a canonical chart. We 
proceed by constructing Pr in exactly the same way as 
Kuchaf did for the Schwarzschild black hole. Then we 
show that the transition to the new chart is indeed a 
canonical transformation by displaying its generator. 

Kuchaf [l3 proposed that (P, P, Pr, Pp) should form 
a canonicalchart whose coordinates are spatial scalars, 
whose momenta are scalar densities and which is such 
that H r (r) generates Diff R. This means that 



H r = t'P t + R'P R 



LP'l 



t'P t + R'P R + F'Pp m 



(4.9) 



Substituting the expressions derived earlier for P and Pf 
into the above one finds 



Pr = Pr - 



LP L LP l 



2R 



2RF RL 2 T 



(4.10) 



where A = (RR')(LP L )' - (RR')'(LP L ). We must now 
show that the transformation 

{t,R,L,P t ,P r ,P l )^(t,R,F,P t ,P r ,P f ) (4.11) 

is a canonical transformation. 
To do this we solve the set 



,JQ](r') , SQ 



poo 

Pt (r)=J2 / dr'Pjir')- , .. , 

^ Jo oqi{r) Sqi{r) 

A JO 



Spi(r) 5pi(r) 



(4.12) 



where (p^ q) and (Pj, Qi) are respectively the old and the 
new phase-space variables and where 0[qi,Pi] generates 



*Pf(t") is the equivalent of Kuchaf 's Pm{t), the momentum 
conjugate to the black hole mass function, M(r). The par- 
allel between our construction for LTB metrics and Kuchaf 's 
construction for the Schwarzschild black hole, which inspired 
this work, is remarkable because the Schwarzschild metric 
(with a varyin g m ass, M{r)) is not diffeomorphic to the 
LTB metric in (2.1) except in the black hole case, i.e., when 
F = 2M = const., and / = 0. 



the transformation. Because we know P and Pp in terms 
of (P, L, P R , Pl), the four non-trivial equations, viz., 



Pl{t) -- 
P R (r) = Pfl(r) + 



dr' Pp(r') 



SG 



8L(r) SL(r) 
dr Pf (r ) - 



d r>P F (r>) 5F ^ 



SR{r) 

sg 



sg 

SR(r) 



SP L (r) SP L (r) 
dr'Pp(r')^l + ^- 



(4.13) 



can be solved for g, and Pr can be recove red u sing the 
second of the above equations. The last of ( 4.13 ) i mplie s 
that g is independent of Pr. The third equation in ( 4.13| ) 
reads 



sg 

SPl 



2P F P L 
R 



LPl 
R 2 



R lZ 



and can be integrated to give 
1 



g = 



dr 



LP L - 

g 2 [R,L] 



Pl 

P 2 



LPr 



(4.14) 



2 1 PP' - LP L 



(4.15) 



and the first equation ensures that ^2 is independent of L. 
We will take_C/2 to be independent also of R and u se g to 
determine Pr fr om t he second equation in (4.13). This 
gives precisely (4.10). Next, we verify that the trans- 
formation has not introduced fresh boundary terms by 
computing the difference between the old and the new 
Liouvillc forms, 



dr[P R SR + P L SL - PrSR - P F 6P] 



dr{8 



LP L - l -RR'ln\ RR, + LPL \ 
2 PP' — LP L 



1 „ , RR' 
2 RSR ^RRT 



LPr 



LP L ' 



(4.16) 



Kuchaf Jl4j has shown that the fall-off conditions (3.7) 
at infinity imply that 



1 -RSRIu\ RR ' + LPl \ 
2 1 PP' - LP L 



LP l SR 

R! 
0(r~ e ) 



(4.17) 



and hence vanishes a t infi nity. Again, as r —* 0, the 
fall-off conditions in ( |3.13|) imply that P -> R (t)r n / 3 , 
SR -> 5R (t)r n / 3 , R' ~> P (t)r"/ 3 - 1 , L -> Laity 1 / 3 - 1 
and P L -> 0{r 2n / 3+1+e ). We find that 



2 l RR> 



LP L 



LPr 



LP L SR 



ra±0 



P' 



.n + l + e 



) (4-18) 



and therefore also v anishes at the origin. The functional 
g defined by (4.15) is well defined. The integrand is of 



G 



order r~( 1+ ^ at infinity and of order r 4n / 3 + 1 + e a t the 
origin, avoiding divergences at both places. 

There are (infinite) boundary terms at the horizon, 
when T = 0. It can be shown, however, that the con- 
tribution from the interior and the exterior cancel each 
other. In the same way, there will be contributions at the 
boundary between the interior of the star and its exte- 
rior or more generally at any frontier between two LTB 
regions described by different mass functions. Again, if 
the mass function is continuous across the boundary and 
regions are consistently matched by equating both the 
first and second fundamental forms, then the contribu- 
tion from one side will cancel the contribution from the 
other. 

In terms of the new variables, the action in ( |3.3[ ) can 
be expressed as 

S = [ dt [ dr[P T T + P R R + P F F- 



NH 9 - N r H 9 r ] + S aE (4.19) 



with 



H = - 



F'T~ X R' + TPfPr 



2L 



H r = t'P t + R'P R + F'P F 



N + (t)M + (t) 



(4.20) 



Let us now re-express the boundary action in a more 
convenient form. As Kuchaf pj] has emphasized, N + (t) 
must be treated as a prescribed function of t. This is to 
avoid the conclusion that the total mass M + as measured 
at infinity is zero, which would follow from varying the 
lapse function at infinity. N + (t) must be chosen and, 
once chosen, held fixed. The freedom in choosing this 
function can be combined with the freedom we have of 
setting the dust proper time at infinity to correspond to 
the parameterization clocks there. The lapse function is 
the rate of change of the proper time with the coordinate 
time at infinity, so it is natural to set N + (t) = f+(t) and 
write the surface action as 



Sdz, 



dt[M+t^ 



(4.21) 



It is linear in the time derivative, f+, and defines a one 
form which can be re-written in terms of the mass func- 
tion, F, as 



M. 



1 Z" 00 

I + 6t+ = --J dr{F5r)' 

1 f°° 

= — / dr[F'5r -t'5F + 8(Ft')]. (4.22) 

2 Jo 



The first two terms on the right hand side may be ab- 
sorbed into the Liouville form of the hypersurface action 
(they modify the canonical momenta) and the last term 
is an exact form which can be dropped. The action is 
thus_expressed entirely as_a hypersurface action. Defin- 
ing P T = P T - F /2 and P F = P F + r'/2, we have 



S= dt 



dr[P T i + P R R + P F F- 
- NH 9 - N r H 9 } 
where the constraints in the new chart read 



(4.23) 



H 



F'T- X B! + T{P F - t'/2)P } 



+ (p r + F'/2)yi + 

t'P t +R'P r + F'P f 



L 2 



(4.24) 



Finally, eliminating the momentum P F from the Hamil- 
tonian constraint by using the momentum constraint, we 
obtain 



(P T + F'/2) 2 + TP R 



l2 



F 



0. 



(4.25) 



This is the simplified constraint referred to in the intro- 
duction. It takes the place of the Hamiltonian constraint 
in ( 4.24 ). In the following section we quantize the con- 
straints and solve the Wheeler-DeWitt equation for arbi- 
trary non-constant mass functions. 



V. QUANTIZATION 

From the dynamical constraint in the previous section, 
one reads off the DeWitt supermetric, jab, on the config- 
uration space X a = (r, R), 



lab 



1 




(5.1) 



To quantize the system, the momenta must be turned 
into operators which act on the state functional. There is 
a standard procedure to do this, which follows a proposal 
due to Vilkovisky ||,|| and DeWitt @: exchange the 
classical momenta for covariant functional derivatives, 



P, 



8X a (r) 



r, 



(5.2) 



where Y is the connection belonging to the configuration 
space metric, 7^. Then defining 

*[r, R, F] - iT F '^Mr)dr^ R n (5 3) 

we see that \& obeys the "Klein-Gordon" equation with 
a potential, 



7 ah V a V b 



pi 2 - 
4F 



SLf[X,F] = 0. 



(5.4) 



The metric in (5.1) is positive definite outside the horizon 
(J- > 0) and indefinite inside (T < 0). The functional 
equation is therefore elliptic outside the horizon and hy- 
perbolic inside. Furthermore, the configuration space is 
flat and this metric is brought to a manifestly flat form 
by the coordinate transformation 



7 



P* = 




(5.5) 



In terms of P*, the Wheeler-DeWitt equation reads 



5 2 S 2 F' 2 
St 2 6R 2 + 4P 



*[r,P,P] =0, 



(5.6) 



or 



5^ ~ F' 



*[r,P,P] = 0, (5.7) 



where the negative sign within the square root refers to 
the region outside the horizon, the positive sign to the in- 
terior and where P* is conjugate to P*. Invariance under 
spatial diffeomorphisms is enforced by the momentum 
constraint, 



J 6 



*SR, 



F'- 



SF 



*[r,P,P] = 0, 



(5.8) 



and, together, eqs. (5J3) and (5^8) define the quantum 
theory whose inner product is given by 

/>OG />00 

(*i,* 2 )= / dR^\^ 2 = ]Tdi?»(r)*I* 2 

JR« (0) JR., (0) r 

(5.9) 

where P(0) represents the physical radius at the center, 
r = 0. The inner product ensures the hermiticity of the 
momentum, P*, conjugate to P*. The norm of a quantum 
state under this scalar product is formally t independent 
provided that h defined in ( |5.7|) is self-adjoint. However, 
this occurs only in the linear sub-space in which the op- 
erator =fP? + F' 2 jAT admits positive eigenvalues. When 
F' = 0, for example for the Schwarzschild black hole, 
this is true of all states with support only in the interior 
of the hole pq| . Such states alone may be ascribed a 
probabilistic interpretation. 

While specific models of stellar collapse must be ana- 
lyzed individually, a solution for general F(r) ^ const, 
can be given as follows. The momentum constraint re- 
quires the wave functional to be a spatial scalar, therefore 
consider a solution of the form 



*[t,P,P] = exp 



drF' [r)W{r{r),R{r),F(r)) 



(5.10) 



where W(r, P, P) is an arbitrary (complex valued) func- 
tion of the coordinates (and not their derivatives) which 
is to be determined. The integrand in the exponent is 
clearly a spatial density because, while r(r), R(r) and 
F(r) are spatial scalars, P'(r) is a density. It follows 
that the wave-functional so d efin ed will obey the mo- 
mentum constraint. Indeed, (|5.8| ) simply requires that 
W admit no explicit dependence on r. It is well-known 
that the Wheeler-DeWitt equation requires some regu- 
larization, involving as it does two functional derivatives 
taken at the same point. We will follow DeWitt |^6| (see 



also p7| ) and require that the delta function vanishes 
in coincidence limit, as does the coincidence limit of its 
derivatives to all orders, i.e., we let 5(0) = = 6^(0) 
in what follows (where the superscipt (n) refers to the 
n th spatial derivativ e of the delta function). Then the 
dynamical equation (|5.6|) reads 



F I2 



dwy fdwy i_ 

dr J \dRj + T 



* = 



(5.11) 



Let us choose W to be of the form W = -ir + U(R, F). 
This ansatz is motivated by the fact that P'/2 actually 
represents the radial mass-ene rgy d ensity within spher- 
ical shells. Inserting it into (5.11) gives the following 
equations for U(R, P) 



dU 



±i 



F 



<9P* "VP-P' 
W = -ir± iVF I dR 



P 



P- P' 



(5.12) 



outside the horizon (P > P), and 



±i 



F 



<9P* 

W = -ir±iy/F I dR 



R 



F-R 



(5.13) 



inside (P < P). In both regions the equations are easily 
integrated. They give 



W = -ir±2iVF 
outside, and 

W = -ir =F 2iVF 



VR- x/Ptanh" 1 \ j — 
R 



VR- v^tanh -1 \l- 
F 



(5.14) 



(5.15) 



inside. 

Any collapse model that approximates reality must in- 
volve two or more regions of differing mass functions (the 
simplest example would be a constant mass function rep- 
resenting the Schwarzschild exterior of a star matched to 
another non-constant mass function representing the in- 
terior of the star). In general the mass function must 
be continuous but it may not necessarily be differen- 
tiable across the boundary between these regions. At 
such boundaries, the wave-functional must be appropri- 
ately matched, i.e., required to be both continuous and 
differentiable. Specific models which exhibit some of the 
interesting classical features whose correct understand- 
ing cannot be had without appealing to quantum gravity, 
as mentioned in the introduction, will be described in a 
forthcoming publication. 

We have shown that there exists a canonical chart 
analogous to that used by Kuchaf to describe the 
Schwarzschild black hole and which describes the 
marginally bound collapse of inhomogeneous, pressure- 
less dust. This chart enjoys several advantages over the 
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original. For one, the proper time of the collapsing dust 
enters naturally and serves as a time variable in the 
Wheeler-DeWitt equation. Again, the chart describes 
the collapse in terms of a more transparent and phys- 
ically meaningful set of variables, yielding a dynamical 
constraint that is greatly simplified and whose solution 
we have displayed over the entire class of models being 
considered. 
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